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5 : Non^canonical scaling of the form oc p^a^, where p. and p. are the energy densities of dark 

I energy and dark matter, respectively, provides a natural way to address the coincidence problem 



-why the two densities are comparable today. This non-canonical scaling is achieved by a suitable 
interaction between the two dark components. We study the observational constraints imposed 



^p^l on such models from the cosmic microwave background (CMB) anisotropy spectrum. Using the 

O 

in ' recent WMAP results for the location of the CMB peaks, we determine the admissible parameter 

■ space and find that interacting models are well consistent with current observational bounds. 
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I. INTRODUCTION 



The picture that emerges from the present cosmological observations portrays a spatially 
lat, low matter density universe, currently undergoing a stage of accelerated expansion 
3]. Within Einstein's relativity, the simplest explanation for this acceleration requires two 
dark components: one in the form of non-luminous dust ("dark matter") with negligible 
pressure, contributing roughly one-third of the total energy density of the Universe and 
clustering gravitationally at small scales, and the other one a smoothly distributed compo- 
nent having a large negative pressure ("dark energy") and contributing about two-thirds of 
the total energy density of the Universe. Although the immediate candidate for this dark 
energy is the vacuum energy (i.e., a cosmological constant A), alternative scenarios, with 
acceleration driven by a dynamical scalar field, called "quintessence", were proposed from 
different perspectives 

One problem that afflicts most of the models proposed so far is the so-called ^^coincidence 
problem" or ^^why now?^ problem j^. The essence of this problem is as follows: as these 
two dark components redshift with expansion at different rates one may ask "w/iy the ratio 
of these two components is of the same order precisely today?, i.e., why = 0{1) f 
There have been different approaches for solving it j^, |5|. The one to be discussed here 
considers some non-canonical scaling of the ratio of dark matter and dark energy with the 
Robert son- Walker scale factor aif) jaBi 

— oc a ^ , (1) 

px 

where pm and are the energy densities of the dark matter and the dark energy, respectively, 
and the scaling parameter ^ is a new quantity which assesses the severity of the problem. 
Thus, = 3 corresponds to the LCDM model and ,^ = to the self similar solution with 
no coincidence problem. Hence, any solution with a scaling parameter < ^ < 3 makes the 
coincidence problem less severe. The standard noninteracting cosmology is characterized 
by the relation ^ = —Zw^, where = Px/Px is the equation of state parameter for the 
dark energy which we here assume to be constant for simplicity. Solutions deviating from 
this relation represent a testable, nonstandard cosmology. Such deviation from standard 
dynamics can be obtained if the two dark components are not separately conserved but 
coupled to each other. This proposal has been explored in 3, IsJ and looks promising as a 
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suitable mutual interaction can make both components redshift coherently. 

Now, it seems advantageous to use the available observational information to constrain the 
nature of dark energy with minimal theoretical input rather than to perform a detailed fit to 

n 

a particular model with a specific potential. This approach has been followed in !(|], where 
high redshift supernovae survey, Alcock-Paczynski test to quasar pairs and evolution of 
cluster abundances have been used to constrain the scaling parameter ^ and the dark energy 
equation of state parameter Wx for the case of separately conserved components. As shown 
in 0] the high redshift SNIa data given in P| cannot discriminate between interacting models 
and the "concordance" ACDM model. Therefore, the target of this work is to constrain the 
parameters ^ and Wx for this kind of non-standard cosmology with the Cosmic Microwave 
Background (CMB) anisotropy spectrum. To this end, we use the Wilkinson Microwave 
Anisotropy Probe (WMAP) [lo| and BOOMERanG data for the location of the peaks 
in the angular spectrum. As it turns out, the scaling interacting cosmology is well consistent 
with the observational bounds. The parameter space corresponding to interacting models, 
seems even to be favoured compared with the parameter space of the standard noninteracting 
models. 



Section 2 succinctly reca! 



the scaling model of Ref. j?]; this component is crucial, of course, if one wishes to test the 
interacting cosmology with the CMB data. In section 3 the positions of the CMB peaks as 
witnessed by WMAP are used to constrain the model. Section 4 summarizes our findings. 
Finally, the Appendix collects the set of formulae employed in our analysis. 

II. SCALING SOLUTIONS 

We consider a Friedmann-Lemaitre-Robertson-Walker (FLRW) universe such that its 
total energy density consists of radiation (whose energy density lies at present by four orders 
of magnitude below that of matter), matter with negligible pressure (that encompasses 
both baryonic and non-baryonic components) and dark energy, 

P = Pr + Pm + Px ■ (2) 



Is and slightly generalizes, by introducing a radiation component. 
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The corresponding pressures are 



Pr = |Pr, Pm « Pm, Px = U^xPx ■ (3) 

The m and x components are supposed to share some couphng so that their energy 
densities obey the balances 

Pr + 4:Hpr = , (4) 
Pm + SHpm = Q , 

p^ + 3H{l + w.^)p.^ = -Q , (5) 

where H = a/ais the Hubble factor and the (non-negative) quantity Q measures the strength 
of the interaction. For simplicity this setup is chosen such that the entire matter compo- 
nent takes part in the coupling. Alternatively, one may treat the baryonic component as 
separately conserved. The constraints obtained below do not depend on whether or not the 
baryons are included in the interaction. 

From Eq. (4) it follows that the radiation redshifts as 

pr = pro (^^) , (6) 

where oq is the current value of the scale factor. 

We are interested in solutions with the following scaling behaviour for the two dark 
components, 

!^ = r,M , (7) 

Px V a / 

where ro denotes the ratio of both components at present time and ^ is a constant parameter 
in the range < < 3. In j^, it was shown that such scaling solutions follow from an 
interaction characterized by 

V - -3ii — — — Pm , (Oj 

1 + ro(l + zj'i 

where 1 + z = ao/a(t). For the standard cosmology without interaction - i.e., when Q = 
-, we have ^ = —3wx- The ACDM model {wx = —1) is the special case with ^ = 3. Any 
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deviation from Q = 0, i.e, (^/3) + Wx = 0, implies an alternative, testable, non-standard 
cosmology. With the interaction (jHJ it is straightforward to find 



n -L n - n lltjl H _L ^^3(1+10,) 
Pm + Px — PmO U + Z) ^ 

Again, for = —I and ^ = 3 this reduces to 

PmO I 



l + ro(l + z) 
1 + ro 



3wx 



Pm + Px = [1 + ^^0(1 + Z) 



ro 



(9) 



(10) 



which is indeed the prediction of the ACDM model. 

With the evolution of the components as given in Eqs.(6) and (9), the Friedmann equation 
can be written as 



8nG 



Hi 



ro 

+ Vlroil + Zf 



1 +ro(l + z)g 
1 + ro 



(11) 



where we have introduced the dimensionless density parameters 



mO 



PmO 
Pc 



rO 



8nG 



We can recast the last equation into 



(12) 



with 



X{z) = (1 + ;2)3(l+-.)-4 



l + ro(l + z)^ 
1 + ro 



rO 



mO 



(13) 



(14) 



The above form of is useful in testing various models against CMB observation jl2i |. 



III. CONSTRAINTS FROM CMB 



The CMB acoustic peaks and troughs arise from oscillations of the primeval plasma 
just before the Universe becomes translucent. These oscillations are the result of a 
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balance between the gravitational interaction and the photon pressure in the tightly bound 
photon-baryon fluid. The locations of the peaks corresponding to different angular momenta 
depend on the acoustic scale Ia_, which in turn is related to the angular diameter distance D 
to the last scattering, and on the sound horizon si^ at last scattering through I a = ttD / Sis 



lo Liie labL bCciLLeiiiig, ciiiu uii Liie bouiiu liuii^on ai 
[isj l . To a good approximation this ratio for I a is 



Ia = n^-^ , (15) 

CsTls 

where f(= / a'^dt) is the conformal time and the subscripts and Is represent the time at 
present and at the last scattering era, respectively, is the average sound speed before last 
scattering deflned by 

~-i f^''" 

Cs = Ti, / Csdf, (16) 







with 



c;' = 3 + !^, (17) 



where pb stands for the energy density of baryons. 

In an ideal photon-baryon fluid model, the analytic relation between the position 1^ of 
the m-th peak and the acoustic scale Ia is = m Ia- But this simplicity gets disturbed by 
the different driving and dissipative effects which induce a shift with respect to the ideal 
position fl3^]. This shift has been accounted for by parametrizing the location of the peaks 
and troughs by 

Im = Uim - (pm) = Uim - - 5(f)m) , (18) 

where is an overall peak shift, identifled with 0i, and = 4>m — is the relative shift of 
the m-th peak. This parametrization can be used to extract information about the matter 
content of the Universe before last scattering. Although it is certainly very difficult to derive 
an analytical relation between cosmological parameters and phase shifts, Doran and Lilley 
[isj l have given certain fitting formulae which makes life much simpler. These formulae do 
not have a prior and crucially depend on cosmological parameters like spectral index (n^), 
baryon density (cjf, = Q^h^), (normalized) Hubble parameter (h), ratio of radiation to matter 
at last scattering (r/s) and also on flf^, representing the dark energy density at the time of 
recombination. We use these formulae (collected in the Appendix) to specify the positions 
of the peaks in the scaling model and constrain its parameter space with the WMAP data. 



It should be mentioned here that although these formulae were obtained for quintessence 
models with an exponential potential, they are expected to be fairly independent of the 
form of the potential and the nature of the late time acceleration mechanism, as shifts 
are practically independent of post recombination physics. It should also be stressed that 
the errors associated with the analytical estimators for the peak positions we are using, 
determined in comparison with the CMBFAST for standard models, is less than 1% jisl ]. 

Let us now turn to discussing the latest available CMB data on the positions of the peaks. 
The bounds on the locations of the first two acoustic peaks and the first trough from the 
WMAP measurement of the CMB temperature angular power spectrum are 

Ip, = 220.1 ±0.8, 

Ip^ = 546 ±10, (19) 
Id, = 411.7 ±3.5; 

notice that all uncertainties are within la and include calibration and beam errors. The 
location for the third peak is given by BOOMERanG measurements jl^ 

Ip, = 8251}° . (20) 

Now, with the help of the above formulae we calculate the position of the peaks in 
the CMB spectrum in scaling models and constrain the model parameters to the values 
consistent with the observational bounds. The acoustic scale is determined in terms of the 
conformal time f in Eq. (15). From Eq.(13) we get 

Tis = dr = -^j^ / -— , 21 

and 

To = / dT= -^j^ / — - , 22 

where X{a) is given by Eq.(14) and we have chosen qq = 1. 

Inserting the above expression in equation (15), we obtain an analytical expression for I a, 



l-A — — 

Co 



/o 



1 da 



X{a) 



rais da - ' (^^) 

JO X{a) 

From the computation of the acoustic scale by Eq. (23), the equations for the peak shifts 
Eq. (18), and the fitting formulae in the Appendix, we look for the combination of the 
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model parameters that are consistent with the observational bounds. We have plotted the 
contours consistent with the bounds of the first three acoustic peaks and the first trough 
corresponding to the WMAP and BOOMERanG data given by Eqs. (19) and (20) in the 
^ — Wx parameter space for different values of and ^^0- The investigated cosmo logical 
parameter space is given by (n^, h,uJh,Qrn)- Throughout this paper, we have neglected the 
contribution from the spatial curvature and massive neutrinos. We have also neglected the 
contributions from gravitational waves in the initial fiuctuations. Because of the rather tight 
WMAP constraint on Ub, {uJb = 0.0224 ± 0.0009 Q), we have assumed uj^ = 0.0224 in our 
calculations. To have a clear idea of the dependence on the parameters we have drawn 
two different types of plots. The first three figures show the — parameter space for a 
particular value of h{= 0.71) and different values of Ug and fimo- The next two figures depict 
the Qmo ~ h parameter space for different values of Ug, and Wx- 

As already mentioned, we have assumed in the general outline of Sec. 2 that the entire 
matter component, including the baryons, takes part in the interaction with the dark energy. 
We have also investigated the case where only the non-baryonic matter part is coupled to 
the dark energy, while the baryonic energy density is locally conserved. The resulting plots 
do not depend on these different ways to implement the interaction. 

In figure 1, we have plotted the contours corresponding to the bounds of the first three 
peaks and the first dip (Eqs. (19) and (20)) in the ^ — parameter space with h = 0.71 
and Qmo = 0.2 for different values of Ug. Figures 2 and 3 represent similar contours with 
same h and Ug but different Qmo (0.3 and 0.4, respectively). To facilitate the analysis, all 
these figures show the non-interaction line, ^ + 3wx = 0. It is apparent that the interaction 
is severely restricted by the CMB data. And simultaneously it is very interesting to notice 
that in Figs 2 and 3 the line ^ + 3^^ = is clearly outside the allowed CMB regions which 
implies that with the assumed priors {h = 0.71 and the mentioned value of Ug and ^mo) 
there is no parameter space at all consistent with a noninteracting cosmology, including the 
ACDM model. 

Since the above important conclusion depends crucially on the priors chosen, we investi- 
gate the same models in a different parameter space. Figure 4 depicts the contours for the 
same peaks and dips in the Qmo ~ h parameter space with Ug = 0.97 and with three different 
values of ^ and Wx- We have chosen the values of ^ as 1.0,2.0,3.0. The higher the value 
of ^, the more acute the coincidence problem. For Wx we have chosen —0.5, —0.75, —1.0. 
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Of course, for ^ = 3.0 and Wx = —1.0 the model collapses to the LCDM model. Figure 5 
represents the same contour as figure 4 with Ug = 1.0 and the same values for ^ and Wx- 

IV. DISCUSSION 

We have investigated the constraints imposed by the observed positions of the peaks of 
the CMB anisotropy spectrum -as witnessed by WMAP and the BOOMERanG experiments 
jl(] |. jll|-. on the non-standard scaling interacting cosmology ^ + 3wx ^ model [3]. The 
most important consequence of our analysis is that the non-standard scaling interacting 
cosmology shows good consistency with the observational bounds for a variety of param- 
eter combinations. The parameter space favoured by the CMB is larger for interacting 
cosmological models than for noninteracting ones. 

In Fig.l, for Q^o = 0.2, the CMB data suggests a constrained parameter space for all 
the four values of Ug, whereas in Fig. 2, for Qmo = 0.3, we get an allowed region only for 
Us < 0.97. In Fig 3, for Q^o = 0.4, we do not have an admissible parameter space at 
all. In Figs. 1 and 2, Wx varies in the limits —0.38 > Wx > —1.08 over the entire range 
oi ^{0 < ^ < 3). Thus, from the whole set of figures it is rather obvious that the scaling 
interacting model favours lower value of and a moderate range of values for Qmo- On the 
other hand, the noninteracting model satisfies the parameter space bounded by CMB data 
only in Fig 1., i.e, for Qmo = 0.2. In figures 2 and 3 the noninteraction line stays well beyond 
the reach of the parameter space allowed by the CMB data. This includes the concordance 
ACDM model as well. Thus, the scaling model is better consistent with the CMB data and 
it is compatible with a larger parameter space than the noninteracting standard model. 

It is worthy of note that for ^ < 1, the region bounded by the contours of Figs. 1, 2 

and 3 reduces practically to a line. For the stationary solution with no coincidence problem, 

^ ~ 0, we practically get a particular value for Wx satisfied by the bounds of WMAP and 

BOOMERanG data. This value of Wx varies only with Q^o -between Wx — —0.38, when 

^mo = 0.2, and Wx — —0.52, when Q^o = 0.4. (It should be borne in mind that claims 

I I 

according to which observations imply that Wx < —0.87 [18^ are based in non-interacting 
cosmologies) . 

As mentioned above in figures 1, 2 and 3, the constrained region crucially depends on h 
and Ug. To make this clear, we present in Figs. 4 and 5 contours in the Qmo — h parameter 
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space for the same ^ and Wx but for different values of ns(0.97, 1.0). From these plots it 
becomes also apparent that for a lower value of ^, i.e., for a less severe coincidence problem, 
we have an admissible f2mo " h parameter space for higher values (lower amounts) of the 
equation of state parameter Wx- For example in the first column of panels in Figs. 4 and 5 
(for ^ = 1.0) we have an allowed parameter space for Wx = —0.50 only, for Wx = —0.75 there 
is no admissible space and for Wx = —1.0 the contours go beyond the range of the parameter 
space chosen. Similarly, as the coincidence problem becomes more acute, i.e., as ^ grows, 
lower values (higher amounts) of become more favoured. This implies that thanks to the 
interaction the dark energy pressure becomes less negative. The bottom-right panel in both 
figures corresponds to the ACDM model (^ = 3.0 and Wx = —1.0). For the panel in Fig. 
4 representing the ACDM model the CMB bounds are satisfied for 0.19 < Qmo < 0.3 and 
h > 0.68, while for the corresponding panel representing the ACDM model in Fig. 5, the 
bounds are satisfied only for h > 0.72 and Qmo < 0.22. This depicts that ACDM favours 
low Us values. This result seems quite consistent with that of Ref . [3| • 



After the recent high redshift SNIa data [la |12| lending additional support to an accel- 
erating universe picture, and with further strings of observations (like those from the SNAP 
satellite) yet to come, one major worry of the community is placed on the nature of dark 



energy 



19|. In this respect, a scaling solution of the type Pm/ Px = ro{ao/a)^ appears to 



be a quite promising tool for analyzing the relationship between the two forms of energy 
dominating our present Universe. And with the available constraints from WMAP and 
BOOMERanG experiments, we can certainly conclude that interacting cosmological models 

+ 3wx 7^ 0) may well compete with the "concordance" ACDM model, they seem even to 
be favoured when compared to the latter one. 

While the present cosmological data are insufficient to discriminate between these models, 
it is to be hoped that future observations of high redshift SNIa as well as other complemen- 
tary data (from galaxy clusters evolution and lensing effects) will decisively help to break 
the degeneracy. 
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APPENDIX 

For completeness, we put together here the formulas used in our search for parameter 
space. These fitting formulae are quoted from the cited literature jis! ]. 

We assume the standard recombination history and define the redshift zis of decoupling 
as the redshift at which th e op tical depth of Thompson scattering is unity. A useful fitting 
formula for zis is given by j20|: 

zi, = 1048 [1 + 0.00124cu,-o-^38][i ^ ^^^92]^ (24) 

where 

0.238ri , on K, ,0.7631-1 „ _nc;cri , 01 i .l.Sli-l 



gi = 0.0783Wfe-"-^^^[l + 39.5cjfe"-^^T > 92 = 0.56[1 + 21.1 



UJt. 



Ub = VLbh^ and uJm = ^m^^- 

The ratio of radiation to matter at last scattering is 

ris = Pr{zis)/Pm{zis) = OMlGu-' {zis/lO') . (25) 

The overall phase shift (f (which is the phase shift of the first peak) is parametrized by 
the formula 

= (1.466 - 0.466n,) [a-^r^^ + 0.29lOfJ , (26) 
where ai and a2 are given by 

ai = 0.286 + 0.626 , (27) 
as = 0.1786 - 6.308 cUb + 174.9 cul - 1168 tul (28) 

The relative shift of the second peak (502) is given by 

S<f2 = Co - cin, - C2rrf' + 0.05 {us - 1), (29) 

with 

Co = -0.1 + (0.213 - 0.123fif,) (30) 
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X exp { - (52 - 63.6QfJ Ub} (31) 

ci = 0.063 exp{-3500 ul} + 0.015 (32) 

C2 = 6 X 10"^ + 0.137 {ujb - Q.Q7 f (33) 

C3 = 0.8 + 2.3nl + (70 - 126ni) Ub. (34) 

For the third peak we have, 

5<^3 = lO-dir^t +0.08(n, -1), (35) 

with 

di = 9.97 + (3.3 - 30fJ Ub (36) 

d2 = 0.0016 - 0.0067Qf, + (0.196 - 0.22QfJ Ub 

^ (2.25 + 2.77QfJxl0-5 ^ ^^^^ 

The relative shift of the first trough is given by 

S(f 3/2 = &o + ^ir/f exp(62n,) + 0.158(n, - 1) (38) 

with 



60 = -0.086 - 0.079Qf^ - (2.22 - 18.1QfJ Ub 

- (140 + 403fif,) Lu'b , 
hi = 0.39 - 0.980f, - (I8.I - 29.20^") ujb 

+440a;6^ (39) 
62 = -0.57 - 3.8 exp(-2365u;ft2) (40) 

The overall shifts for the second and the third peaks and for the first trough are + 502, 
+ 503 and + 503/2, respectively. In the above expressions, Qf^ is the average fraction of 
dark energy before last scattering, which is negligibly small in the cases discussed here. 
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FIG. 1: Contour Plots of the first three Doppler peaks and the first trough location in the {i,Wx) 
plane with = 0.2 and h = 0.71 for different values of Ug. Black, red, blue and green lines 
correspond to the observational bounds on the first, second, third peaks and the first trough, 
respectively (given in equations (20) and (21)). The upper line corresponds to the non-interacting 
case {^ + 'iwx = 0). To avoid confusion, we plot only the contours relevant for the admissible range. 
For each panel contours not shown are consistent with this range. 
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FIG. 2: Same as Figure 1 with h = 0.71 but r^^o = 0.3. 
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FIG. 4: Contour Plots of the first three Doppler peaks and the first trough location in the (rimO; h) 
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plane with = 0.97 and for different values of and Wx- Black, red, blue and green lines correspond 
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